Strong coupling superconductivity due to massless boson exchange 
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We solve the problem of fermionic pairing mediated by a massless boson in the limit of large 
coupling constant. At weak coupling, the transition temperature is exponentially small and super- 
conductivity is robust against phase fluctuation. In the strong coupling limit, the pair formation 
occurs at a temperature of the order of the Fermi energy, however, the actual transition temperature 
is much smaller due to phase and amplitude fluctuations of the pairing gap. Our model calculations 
describe superconductivity due to color magnetic interactions in quark matter and in systems close 
to a ferromagnetic quantum critical point with Ising symmetry. Our strong coupling results are, 
however, more general and can be applied to other systems as well, including the antiferromagnetic 
exchange in 2D used for description of the cuprates. 
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I. INTRODUCTION 

Strong coupling superconductivity due to the in- 
teraction between electrons and lattice vibrations has 
been successfully studied using the coupled Eliash- 
berg equations! for the frequency dependent normal and 
anomalous self energies of superconductors 2 * 3 *-! The the- 
ory finds its justification in the weakness of the cor- 
rections to the electron-phonon vertex, caused by the 
small ratio of the electron and ion masses*!. This the- 
ory inspired numerous efforts to describe superconduc- 
tivity caused by other bosonaViSiWWWV 6 , even 
though its justification turns out to be considerably more 
subtle in some of those casesi 7 **^**. Important progress 
has been made in the study of pairing due to the ex- 
change of bosons that are collective excitations of the 
fermiono 6 i 7 i 8 iV i 11 i 12 i 13 i 14 i 15 i 16 . In this context, the inter- 
play between superconductivity and quantum criticality 
is particularly interestingSSiSLSSiSiSii 2 ^ as superconduc- 
tivity in correlated electron systems often occurs in the 
proximity of a quantum critical point (QCP) 27 i 28 i 29 i 30 i 31 . 
At a QCP, the pairing boson becomes massless, and 
new and unexpected behavior emerges 2 ^**^**^. A related 
problem occurs in the theory of quantum chromodynam- 
ics at high density where single-gluon exchange becomes 
dominant* 3 * 3 .. The exchange of gluons is believed to cause 
color superconductivity 3 ^ 3 ^* 3 ^. As was pointed out by 
Son-£, and later in Refs. 37-3 — , the color magnetic interac- 
tion in high density QCD is unscreened at low temper- 
atures, i.e. the pairing is mediated by a gapless boson. 
The pairing problem then becomes formally very similar 
to superconductivity at a QCP, even though the transi- 
tion temperatures may be different by a factor as big as 
10 12 . 

In previous studies of the pairing problem near a QCP, 
the authors of RefsiSS* 2 ^* 3 -^ assumed that the effective, 
boson-mediated fermion-fermion interaction u is much 
smaller than the fermionic bandwidth, W (which is gen- 
erally of the same order as Ep), i.e., 
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The limit g -C 1 is often called weak coupling. This 
notation is not quite correct, as near a QCP the smallness 
of g doesn't imply that the system behaves as a weakly 
coupled Fermi liquid - the mass renormalization due to 
the exchange of a gapless boson is still singular in D < 
3 and destroys the Fermi liquid behavior at the QCP 
(see below). To simplify the notations, we nevertheless 
refer to g -C 1 as weak coupling and g » 1 as strong 
coupling. With this notation the Eliashberg theory for 
superconductivity due to electron-phonon interaction* 2 * 3 *-* 
is in the "weak coupling" limit since g oc ^A is small. 
This is due to the smallness of the ratio — of the sound 
and Fermi velocities, while the product A = ppV cp of the 
electron-phonon interaction V ep and the density of states 
at the Fermi level pp can be of order unity. 

The condition j < 1 implies that the pairing comes 
only from fermions in a tiny range of momenta around 
the Fermi surface, i.e., that the system behavior at en- 
ergies comparable to Ep is irrelevant for the pairing. 
This makes the pairing problem universal and allows 
one to use well-established computational techniques, e.g. 
the Eliashberg theoryi. However, for various systems 
of interest the interaction is not necessarily small. In 
particular, the same interaction that leads to the pair- 
ing is often also responsible for the onset of order at 
the QCP. Generic density-wave instabilities come from 
fermions with energies O(Ep) and require g to be of or- 
der unity 3 ^. In the cuprate superconductors, to which 
the ideas of collective-mode mediated d— wave pairing 
was appliecUi^iS* 2 ^, the Hubbard interaction U is at least 
comparable to W as is evidenced by e.g. the Heisenberg 
antiferromagnetism at half-filling. For color supercon- 
ductivity, the effective coupling u is also not necessarily 
small compared to Ep, and g well may be larger than 1. 

These arguments call for an understanding of the pair- 
ing problem beyond the "weak coupling" limit. In the 
present paper we extend previous "weak-coupling" stud- 
ies of the pairing mediated by a gapless boson to the truly 
strong coupling limit j » 1. For definiteness we con- 
sider pairing of 3D electrons mediated by a scalar boson 
which is gapless at q = 0. This model describes p-wave 
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FIG. 1: Schematic phase diagram for the superconducting 
transition temperature T c and the pairing instability temper- 
ature Tpair as function of the dimensionless coupling constant 
g. While T c ~ T pa i r for weak coupling, an intermediate regime 
T c < T < Tp a i r with phase (and amplitude) fluctuations oc- 
curs in the strong coupling limit. 

superconductivity near a ferromagnetic Ising QCP, and 
color superconductivity of quarks. However, the results 
at strong coupling are quite general and can be applied 
to other systems as well, including the antiferromagnetic 
exchange in 2D used for description of the cuprates. We 
discuss applications to other systems in a separate sec- 
tion. There is also a connection between our model and 
the interaction between conduction electrons, mediated 
by transverse photons ^2*^. However, as shown in&, the 
exchange of transverse photons does not lead to super- 
conductivity. 

A word of caution. In the context of quark pairing 
mediated by gluons, the equation for the pairing vertex 
has only been derived in a gauge invariant manner at 
weak coupling, — . At strong coupling, antiparticle pair- 
ing, neglected in our model, may come into play. Still, 
qualitatively, the results obtained assuming only parti- 
cle pairing likely remain valid at both weak and strong 
coupling. 

The main results of this paper are summarized in the 
phase diagram, Fig[I] In the weak coupling limit, we 
find, in agreement with Son^ and others^L^, that the 
transition temperature behaves, to leading exponential 
order, as 

log— = — , 2 

where loq ~ Ep/g. This result is parametrically larger 
than the usual BCS result: logwo/T c oc 1/g, and the 
difference is due to the gapless nature of the pairing boson 
(see below). Still, T c is exponentially small at small g. 

At g = 0(1), T c becomes of order Ep, although the 
prefactor for T c is a small number. At even larger g, we 
find two characteristic temperatures. The larger temper- 
ature, Tp a i r , sets the onset of pairing, and is of order Ep 



(again, with a small prefactor). The smaller temperature, 
T c is of order y/ujoEp ~ Ep/^/g <c Ep. This temper- 
ature is determined by the superfluid stiffness, and sets 
the scale for phase coherence, i.e., of the actual supercon- 
ductivity. In between T c ~ Ep/ \fg and Tp air ~ Ep, the 
system displays pseudogap behavior: pairs of fermions 
are already formed, but do not move coherently. 

The accuracy of the computations is a central is- 
sue for the theoretical analysis near a QCP. At small 
g, one can use Eliashberg theory since the relevant 
bosonic and fermionic frequencies are much smaller than 
Ep. Although the frequency-dependent self-energy is 
not small due to the near-criticality, vertex corrections 
and the momentum-dependent self-energy are exponen- 
tially small (see below). At g = 0(1), typical frequen- 
cies become of order Ep, vertex corrections become 0(1), 
and the momentum-dependent self-energy becomes of the 
same order as the frequency-dependent self-energy. In 
this situation, no reliable theoretical scheme is possible. 

A naive expectation would be that at strong coupling, 
vertex corrections get even stronger. We show, however, 
that at larger g, vertex corrections actually saturate at a 
value 0(1) and do not grow with g. At the same time, the 
momentum dependent term in the self-energy again be- 
comes small compared to its frequency dependence, this 
time the relative smallness is in 1/g. Furthermore, the 
pairing problem at g 3> 1 still involves fermions with en- 
ergies below Ep for which the density of states can be 
approximated by a constant. As a result, the new version 
of the Eliashberg theory (more accurately, the local the- 
ory) becomes qualitatively valid. This new local theory is 
different from the original Eliashberg theory in that the 
lattice cannot be neglected, and the existence of a finite 
bosonic bandwidth now plays a crucial role. Still, like in 
the Eliashberg theory, we derive closed-form equations 
for the fermionic self-energy and the pairing vertex. To 
justify the local approximation at g > 1 quantitatively, 
we extend the model to N fermion flavors and consider 
the limit of large N— . In this case, vertex corrections 
become of order 1/N and can be safely neglected. 

The structure of the paper is as follows. In the next 
section we set up the model and define the large N limit 
used to perform the strong coupling calculation. In Sec. 

III we briefly discuss the weak coupling limit and present 
an alternative derivation of Son's result 32 for T c . In Sec. 

IV we solve the pairing problem at strong coupling. In 
Sec. V we analyze gap fluctuations and demonstrate the 
existence of two characteristic temperature scales. In Sec 
VI we justify our computational procedure. In Sec VII 
we discuss other systems, including the cuprates. The 
last section presents our conclusions. Several technical 
details are presented in the Appendix. 



II. MODEL AND LARGE N EXPANSION 

We consider the pairing problem in which 3D fermions, 
ipkj interact via exchanging a massless, bosonic mode 
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with a static propagator D, 
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Here u> is the effective interaction (with the dimen- 
sion of energy), and kp is the Fermi momentum. Note 
the overall sign in © is opposite to that in systems with 
Coulomb interaction 2 -. Another energy scale in the prob- 
lem is the fermionic bandwidth W (roughly, the scale up 
to which the fermionic dispersion 6k can be linearized 
around Fermi surface, Ek = vp(k — kp)). The ratio of 
the two characteristic energies defines the dimensionless 
coupling constant g in Eq. (Q. 

The interaction © leads to pairing, and also gives 
rise to fermionic and bosonic self-energies, Efc(w) and 
Ilq(O), respectively. The two self-energies are related to 
fermionic and bosonic propagators via 



G fc 1 (w) = iuj — vp(k — kp) — Sfc(w) 

D- X (n) = q 2 + ii q (n). 



(4) 



In order to perform a controlled calculation at strong 
coupling, we generalize the model of Eq. (0} to N fermion 
flavors and rescale vp — > vpN and u —w,N. In what 
follows, we assume that the new vp and u are constants, 
independent on Nm^ 

The generalized Eliashberg theory is a set of three cou- 
pled integral equations for the pairing vertex $ and the 
self-energies E and II. We will primarily be interested 
in the onset of the pairing and consider the linearized 
equation for and normal state expressions for E and 
II. Then the system of three coupled equations has the 
form 



$(lj) = 

EH = 
U q (Q) = 



D q {lu - lu') $ {(J) 



xG kF+q (u')G kF+q (-u'), 

^ / D q {u-uj')G kF+q {u'), 

kF J q ,u' 

^- / G k (uj)G k+q {uj + Q). 

k F Jk.u 



(5) 



We used the notation f ... = f tI^TVV ... with 
Matsubara frequencies Q n — 2wkT and Lu n — 
(2n +1) 7rT for bosons and fermions, respectively. 

Eqs. (JSJ) neglect vertex corrections and the momen- 
tum dependence of E and $. We will argue below that 
at large N, both approximations hold both at weak and 
at strong coupling. Physically, these approximations are 
based on the (verifiable) assumption that bosons are slow 
modes compared to fermions. This allows one to factor- 
ize the momentum integration in Eqs.©- Namely, for 
every given along the Fermi surface, the integration 
over the component q± transverse to Fermi surface in the 
equations for $ and E involves only fast fermions, while 
integrating over the remaining two momentum compo- 
nents gy in the bosonic propagator, one can set q± = 0. 



This implies that the boson propagator actually only ap- 
pears in Eqs. JSJ through the "local" interaction 



d(Q) 



q»dq\\D c 



=o(n), 



(6) 



where go ~ kp is the upper cutoff in the integral over gy. 
As a result, the equations for $ and E in 10) reduce to 

= 3_g I ,jp^) 

E(lj) = -iy J dOsign (0 + w) d (Cl) , (7) 
where the factor of | is for further convenience and the 



coupling constant g is given as 

u 



24ir 2 E F ' 



(8) 



with E* F = vpkp/2. Without further approximation we 
can explicitly solve for II (f2) and find24 



where 



n( g ,n) = 7 H, 



7 = 12-K 2 kl —. 
1 Ep 



(9) 



(10) 



Note that 7 does not depend on N, despite the fact that 
the Landau damping term contains a flavor index N as 
an overall factor. The iV-independence of 7 is the result 
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of our rescaling: in rescaled variables 7 — * N 
finite. 

For the "local" interaction @), we obtain from Eqs. 
@ and ©: 



(11) 



where luq is the characteristic frequency of the bosonic 
degrees of freedom: 



LUq 



Q 3 o 



and we introduced 



Ep 
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12 \irkp 



(12) 



(13) 



Below we will refer to Ep as to Fermi energy. We should 
keep in mind however that our Ep depends on the choice 
of the upper momentum cut off go and is only of the 
same order of magnitude as a actual Fermi energy of the 
system. 

Substituting Ea. l|llJ) into Eq.© and integrating over 
frequency, we obtain 



iT, (lu) = Log [—r log + log — n — 



\lu\ < LUq 



sign (u) guj log ^ \u\ > lu 



(14) 
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We see that luq sets the scale at which the momen- 
tum cutoff in the bosonic propagator begins affecting the 
fermionic self-energy. At low energies, the cutoff is ir- 
relevant, and the self-energy has the form typical for a 
marginal Fermi liquid 45 . At u> > luq, the self energy al- 
most saturates and only logarithmically depends on fre- 
quency. At weak coupling, luq = Ep/g > Ep, and the 
crossover is meaningless as Eq. (|14H only holds up to 
lu r~j Ep (we recall that in obtaining Ea. (|f 4[) we approxi- 
mated the density of states by a constant). The marginal 
Fermi liquid behavior then extends all the way up to Ep. 
At strong coupling luq <C Ep, and the crossover in 
occurs well below Ep. In this situation, marginal Fermi 
liquid behavior only holds at small frequencies ui < Ep/g, 
while at Ef/q < lu < Ef, ~S(lu) depends logarithmically 
on frequency (see Fig. ITTf) . 

The crossover in the self-energy at strong coupling par- 
allels the crossover in the "local" bosonic propagator d(u>) 
in Ea. lfrU) 



d(Q) = 



I log ^ M«"o 
5 ft 



(15) 



Like for the self-energy, this crossover is meaningful only 
at strong coupling, when luq < Ep. 

Substituting the self-energy and d(Q) into the equation 
for the pairing vertex, we obtain 



$(w) 



dJ- 



$(w') 



•i£(u/)| 



log 1 



u 



Strictly speaking, we have to evaluate this equation at 
finite T, because the linearized equation for $ is only 
valid at the onset temperature for the pairing. By rea- 
sons that we outline below, we label this temperature as 
Tpair rather than T c . As we will only be interested in the 
order of magnitude estimate for T pa i r , we adopt a sim- 
plified approach, and instead of performing the discrete 
Matsubara sum, use Ea. lfTBI) at finite T, but introduce 
a lower frequency cutoff at u> ~ T. In the weak cou- 
pling limit, this procedure was shown earlier— to yield 
the same T pa i r (modulo a numerical prefactor), as one 
would obtain by performing an explicit summation over 
discrete Matsubara frequencies. In Appendix B we show 
that the same holds for large g. With this simplification 
we have to solve 



$(w) 



'T palr 

with bosonic kernel 



00 $ (lu 1 ) 

^wM^)\ K{ ^ h (l7) 



1 



log 



CU 



( 18 ) 

In what follows we solve this equation, first in the weak 
coupling limit g -C 1, where we reproduce the results of 
Ref— , and then in the strong coupling limit g> 1. 
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FIG. 2: Characteristic energy scales for the weak (g <C 1), 
intermediate (g ~ 1) and strong (g S> 1) coupling limit. 
While at weak coupling, luq ~ Epjg is large compared to Ef 
and thus irrelevant, it emerges as a new low energy scale in 
the strong coupling limit. 



III. PAIRING PROBLEM AT WEAK 
COUPLING 

At weak coupling, one obviously expects T pa i r to be 
much smaller than luq (see Fig[nJ|. This in turn implies 
that only frequencies lu <C luq are relevant. For these 
frequencies, the self-energy £(w) and the kernel K (lu, lu') 
in 118|) can be simplified to 



£ (lu) = —iLug log 
K(lu,lu') = log 

Eq. (|17J) then becomes 



M 

UJ 



$ (lu) = g 



$ (lu') log 



duj'- 



2 i 



u'(l+gu' log J) 



(19) 



(20) 



This equation yields T pa i r for the pairing in a marginal 
Fermi liquid. Eq. (|2f)|> was solved numerically in Rcf. 44 . 
We show that an analytic solution is also possible. Our 
computational procedure is similar to the one used by 
Son2£. In addition to the approach of Rcf.— we also 
analyze the pairing susceptibility. 

If the two logarithmic terms in the r.h.s. of Eq. H2U|) 
were absent, the equation for the pairing vertex would 
be the same as in BCS theory 4 ^, and Tpair would scale as 
cu e 



-i/g 



However, as Son demonstrated—, the presence 
of the logarithm in the pairing kernel substantially en- 
hances T pa i r at weak coupling and changes its functional 
form to T P air oc LUQe^ 7 *^ 2 ^. The easiest way to see 
this is to introduce logarithmic variables: x = log(^), 
x' = log (^), x T = log (^), and re- write Eq. (|5UJ| with 
logarithmic accuracy as 

,x'$(x') 



$(x) = g 



dx 



1 + gx' 



gx 



XT dx'^l 



gx' 



(21) 



Differentiating both sides of Eq. 1]21H over x, we find 



d$ (x) 
dx 



$(x) 
1 + gx ' 



(22) 



Differentiating one more time, we find that the integral 
equation for the anomalous vertex reduces to a second- 
order differential equation: 



dx 2 



= -.9 



(23) 



The gx term in the r.h.s. of Eq. \\'2'M is due to the 
fcrmionic self energy. We assume and verify afterwards 
that gx -C 1 for all relevant x, and drop this term from 
11231 . The solution of Eq. i|23|) is then elementary: 

$ (x) = Acoa(y/gx) + Bsin(y/gx) . (24) 

The two boundary conditions 



$ (x 
d§(x 



0) 



dx 



0. 





(25) 
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FIG. 3: Particle-particle response function \vv ( w > T) f° r 9 = 
0.1 as function of T for various energies. Insert: \vv ( w > T) as 
as function of energy for T = 1.001T C . 



follow from (t2"T)l and J22J), respectively. They yield A = 0, 
and 



cos^xt) = 0. 



(26) 



The onset temperature T pa j r corresponds to the smallest 
Xt that satisfies (|2*B)l . i.e. to xx^/g = tt/2. Ignoring pre- 
exponential factors, we then reproduce Son's result: 



(27) 



The relevant value of x are xt oc 1/^/g, i.e., <?a; ~ 
gxT ~ ^/ff *C 1. This justifies dropping the <?a; term 
(i.e., fermionic self-energy) from l|23|) . The first order 
correction to Tpair due to the self-energy was analyzed in 
Ref. 46 in the context of color superconductivity. 

To get more insight into the pairing instability, it 
is also instructive to analyze the pairing susceptibility 
Xpp(u;, T). This is done by adding an infinitesimally small 
external pairing field $0 to the r.h.s. of H2Q(! 47 . The pair- 
ing susceptibility is 



X PP (w,T) 



d®(u},T) 



$(w,T) 



At T > Tp a i r , $(w) oc $0; and the pairing susceptibil- 
ity is finite. If the transition is of second order, Xpp( w i T) 
diverges at T pair . In BCS theory, $(w,T) does not de- 
pend on frequency, and x PP (u,T) = glog(T 1 /Tpair) ■ This 
pairing susceptibility is obviously positive above T pa i r , 
diverges at T pa i r for all u, and is negative below T pa ; r , 
implying that the normal state is unstable against pair- 
ing. In our case \vv can easii y be obtained from the 
solution of ilL'31) by changing the boundary condition at 
x = to $ [x = 0) = $0 . We then obtain $(x) given by 
Pljl with 



and 



X PP (uJ,T) 



cos (y/g log |) 
cos {y/g log ^) 



(29) 



y4 = $0) -B = *o tan( 1 /ga;T) 



(28) 



Note that T and wo are lower and upper limits of the 
integration over uj in H2Q(I , hence the pairing susceptibility 
is only defined in the interval T < lo < u)q. We see from 
(|29|l that at the upper boundary, to — u>o, Xpp = 1 at any 
T. This is a clear distinction to the BCS limit. As long 
as T > T pa ; r , the pairing susceptibility remains positive 
everywhere in the interval T < lo < lu despite the fact 
that the solution of the differential equation 1231) for $ (lo) 
is formally an oscillating function of frequency. At T pair , 
log Tp- — aTTf ) anc ^ Xpp diverges for all lo, except to = loq. 
Below T pa i r , Xpp is negative at low frequencies, implying 
that the system is unstable towards pairing. We show 
the behavior of x PP (a ;, T) as function of temperature and 
frequency in Fig lllll 



IV. PAIRING PROBLEM AT STRONG 
COUPLING 

We next analyze the strong coupling limit j » 1. In 
distinction to the weak coupling regime, we now have 
two characteristic energy scales in the problem, loq — 
E>f/9 <C E f , and E F , which is the ultimate upper cutoff 
in the theory (see Fig[nJ|. The issue then is which of the 
two scales determines the onset of the pairing. 

Suppose momentarily that only frequencies lo < loq 
contribute to the pairing. At a; < u>q, the pairing ker- 
nel and the self-energy can still be approximated by Eq. 
{JTJJ, and the equation for the pairing vertex can still be 
reduced to the differential equation (|23|l . In distinction 
to the weak coupling case, however, the term gx, coming 
from the self-energy, is now the dominant term in the de- 
nominator in the r.h.s. of l|23|l . Leaving only this term, 
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we arrive at 



d 2 $ (x) _ $ (x) 
dx 2 x 



(30) 



Note that g drops from this equation because of cancel- 
lation between g factors in the effective interaction and 
the self-energy. 

The solution of Ea.QZfy with $ (x = 0) = is $ (x) oc 
\fxJ\ (2y/x), where J\ is a Bessel function. Substituting 
this solution back into (|21() and assuming that the upper 
limit in the frequency integral in l|21l) is still xt (i.e., 
that only uj < uj are relevant for the pairing) , we obtain 
xt = 3.670(5). This leads to T pair ~ 0.025wo, i-e., to a 
pairing instability at a temperature which is a fraction of 

UJQ. 

This result is similar to McMillan's T pa i r ~ 
u}p,e.~^ 1+3 ^ 9 ~ ^>d for strongly coupled phonon super- 
conductors^ (uip, is Dcbye frequency). However, like for 
phonons, there is actually no reason to restrict the fre- 
quency integral to u> < loq ~ Epjg, since for strong 
coupling there also exists a wide frequency range ujq < 
uj < Ep where, on the one hand, the pairing kernel and 
the self-energy are different from (|19|l . and, on the other 
hand, typical frequencies are still below Ep, i.e., a low- 
energy description is at least qualitatively valid. The 
existence of this extra range raises the possibility that 
the onset of pairing may occur at a temperature of or- 
der Ep, not of order ujq ~ Ep/g. Note that for the 
electron-phonon case, the scale which sets the ultimate 
upper cutoff for the pairing (the analog of Ep in our case) 
is u D Jg&MSW 

To verify whether T pa i r scales as Ep, not as uj , we an- 
alyze the equation for $(cj) assuming that all character- 
istic frequencies are larger than ujq. At these frequencies, 
the pairing kernel and the self-energy are given by 



K(uj,uj') 



-zsign (uj) guj log — 
wo 



2 



1 



1 



\uj-uj'\ \lu + lu'\ 



(31) 



Now the pairing kernel scales as 1/uj, while the self- 
energy is nearly a constant, and only logarithmically de- 
pends on frequency. 

Substituting the pairing kernel and the self-energy into 
(|17|l and using the fact that for all uj < Ep the self-energy 
£(w) exceeds the bare uj, we obtain 



duj'§ (uj') 
21og^ 



\uj-uj'\ \uj + uj'\ 



(32) 



The logarithmic divergence in the r.h.s. of at uj = UJ 
can easily be regularized as the 1/\uj — uj'\ form of the 
kernel is only valid at \uj — u>'\ > ujq. 

We see that the dimensionless ratio T/Ep is the only 
parameter in Ea. (|32|l . except for the log — — term in the 
denominator in Ij32(l . Hence, if this equation has a solu- 
tion at some finite value of this parameter, the pairing 
instability should occur at T ~ Ep. 



The analysis of Eq. (|32|l requires special care because 
of the interplay between the 1/uj dependence of the pair- 
ing kernel and logarithmic behavior of the self-energy. 
The discussion is somewhat technical, and we moved it 
into Appendix A. We find there that the solution of i':>2l) 
at frequencies between T and Ef is 



$(w) 



A—= cos 



/31og 



E F 



(33) 



where (3 = 0.7923(2) is determined from the solution of 
a transcendental equation, and A, <fi are real constants. 

Like at weak coupling, the two limits of the integration 
over uj' in Eq. i1-!2l) for $ imply two boundary conditions 
for Q(uj) from l|33|l . One of them determines the phase 
<p, while the other determines the pairing instability tem- 
perature (the overall factor A in (|33|l cannot be deter- 
mined from the linearized gap equation). For a simple 
estimate of T pa ; r , we use the same boundary conditions as 
in the weak coupling limit, i.e., (i) assume that frequen- 
cies larger than Ep are irrelevant for the pairing and set 
= Ep) = 0, and (ii) assume that d U=T pa)t = 0. 
We then obtain <f> = -| and 

T pair = E F e~^ ~ 0.0676i; F , (34) 



where 1//3* = (l//3)arccos 



^TT^J ~ o^T' In Ap " 
pendix B we demonstrate that the same result, modulo a 
numerical prefactor, is obtained by solving explicitly the 
linearized gap equation for discrete Matsubara frequen- 
cies. 

We see therefore that at strong coupling, the pair- 
ing instability temperature r pa ; r is indeed of the order 
of the Fermi energy Ep, although numerically it is still 
much smaller than Ef. This temperature is larger by 
a factor g than the McMillan-type estimate, T pa i r ~ ujq, 
which ignores the pairing interaction at energies larger 
than the characteristic bosonic frequency. We empha- 
size that in order to obtain Eq. (|34f> . it was crucial that 
we included into consideration the normal state self en- 
ergy renormalization. Had we ignored it, an oscillating 
solution for $(uj) at uj > ujq would not have been pos- 
sible, i.e., no pairing instability would occur at T > ujq 
(see Ref*^). Alternatively speaking, T pa i r ~ Ep is the 
result of the interplay between a non- Fermi liquid behav- 
ior of the fcrmions caused by the logarithmic self energy 
£ (uj) oc gcJo log j^- , and a retarded pairing interaction 
governed by a "local" boson susceptibility d(w) oc — . 

Since Tp a i r ~ Ep, it is inevitable that the magnitude of 
r pa i r is affected by the system behavior at high energies, 
i.e. at lattice scales in the condensed matter context. We 
assumed above that the fermionic density of states is a 
constant. This is indeed only approximately valid at u> ~ 
Ep. To determine T pa i r beyond an order of magnitude 
estimate, one then needs to solve the full microscopic 
problem. Still, lattice effects only modify the prefactor 
in T pair ; the relation T pair ~ Ep is generic and survives 
lattice corrections. 
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V. THE ROLE OF GAP FLUCTUATIONS 
A. Phase fluctuations 

In the weak coupling limit it is known that the tran- 
sition temperature, determined from the linearized gap 
equation, coincides with the temperature where global 
phase coherency sets in. This can easily be seen by eval- 
uating the phase stiffness p s defined as 



E, 



phase 



p s / d 3 x (V</>) 



(35) 



At weak coupling, p s ~ Epkp. Eq. I|35[l can then be con- 
sidered as the continuum limit of an XY- spin model on 
a three dimensional lattice with lattice constant ~ k^ 1 
and exchange interaction ~ Ep. Fluctuation effects in 
this model become effective at temperatures comparable 
to the exchange interaction, i.e., at temperatures com- 
parable to the Fermi energy. Since T pa j r <C Ep, fluctu- 
ations at T ~ Tpai r are ineffective, and phase coherency 
is established as soon as Cooper pairs are formed, i.e., 
T ■ — T 

x pair — x c- 

Consider next the strong coupling limit, where T pa i r ~ 
Ep. In what follows we argue that at strong coupling, 
Ps/kF ~ Ep/ \fg <C Tp a i r . In this situation, phase fluctu- 
ations become relevant well below the onset of the pair- 
ing, and by conventional reasoning52iS&§i ; phase coher- 
ence sets in at 



T c ~ Ep j y/g <C Tp a i r ~ Ep. 



(36) 



This new energy scale is the characteristic energy of a 
boson in the gaped state below T P air- In between T pa i r 
and T c , the systems displays a pseudogap behavior: the 
density of states develops a maximum at a finite fre- 
quency (the tunneling gap), and the spectral weight is 
transformed from frequencies below the gap to frequen- 
cies above the gap. However, the superconducting order 
parameter only develops at T c . 

We now show how we arrived at p s /kp ~ Ep j \fg. The 
superfluid stiffness at T = is obtained by evaluating the 
sum of fcrmionic bubbles made of normal and anomalous 
Green's functions, and is given by 



Ps 



duj- 



$ 2 (w) 



(ujZ (uj)f + $ 2 (w) 



3/2 ' 



(37) 



where p a s ~ Epkp is the stiffness of a BCS supercon- 
ductor. <I>(w) is the pairing vertex at T = and we 
introduced 



(38) 



Using the relation between $(w) and the gap function 
A (uj) = i °ne can write Eq. (|57l) as 

roc duo A 2 (oj) 



Ps = Ps 



o Z(u) [uj 2 + A 2 (w)] 3//2 



(39) 



For a BCS superconductor, Z = 1, and A does not de- 
pend on frequency. The frequency integration in l|39(l 
then yields p s = p®, independent on A. This essentially 
implies that at T = 0, the superfluid density equals the 
full density. 

To obtain p s at strong coupling, we need to know 
A(uj,T = 0) and Z(u,T = 0). The gap A(uj,T = 0) = 
A(w) is obtained by solving the nonlinear gap equation 



AM 



2 



(a(uj') - A(uj)^) d sc (uj - uj') 
^Juj' 2 + A{uj'f 



dJ. 



(40) 

where d sc (fl) is the "local" boson propagator in a su- 
perconductor. In the normal state, d(Q) is given by Eqs. 
(|f III5fl . In the presence of A, the bosonic spectrum itself 
changes due to feedback from the gap opening, and the 
Landau damping transforms into II (f2) ~ 
(Ref^i). This leads to 



7 _0l 

I qE F 



(fi) = 3 log 



0,sc 



n 2 



where 



bJQ,sc 



V^oEp ~ Ep/y/g 



(41) 



(42) 



is the characteristic energy of the bosons in a state where 
fermions are gaped - it has the same physical meaning 
as u>o above T pa i r - For frequencies Ep > > wp,sc we 
have 



d sc (fl) 



w „ 



(43) 



Despite the 1/ft 2 dependence of d sc {tt), the integral over 
uj' in Eq. I|4(J|) remains convergent since the numerator 
vanishes at uj — uj' . We can then safely use Eq. 
for d sc (Q) and drop the restriction that this form is only 
valid above ujq.sc- The gap equation then contains only 
Ep as the energy scale. Accordingly, A(uj) can only be of 
order Ep, if the gap equation indeed has a solution. We 
verified that the solution of (|43|l does indeed exist and 
yields A(uj) = E F f(oj/E F ). 

The expression for Z(uj) follows from the formula for 
the self-energy 



d (uj — uj') uj'duj' 
^Juj' 2 + A{iJf' 



(44) 



Here the restriction that Eq. I|43(l is only valid at frequen- 
cies above ujQ tSC becomes crucial, otherwise the integral 
over uj' in Eq. I|44|l would diverge. Beyond this, the eval- 
uation of the integral is straightforward, and we obtain 



Z(uj < Ep) 



,V2 



, Z(uj^Ef)^1. 



Substituting this Z into (|39|) . we find 

p s {T = 0) ~ p° s /y/g ~wo, S cfef 



(45) 



(46) 
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We see that p s (T — 0)/kp is much smaller than T pa ; r . 
The exchange constant of the XY- model (|35|l is therefore 
^o,sc "C Ep. This leads to our estimate of T c in Eq. (|36(l . 
This estimate is further supported by the fact that a finite 
T, we found that the leading temperature dependence of 
the stiffness varies as a function of 



Wo, a 



i.e., thermal 

corrections to the stiffness indeed become relevant at T 



B. A relation to Eliashberg theory 

We emphasized above that our strong coupling theory 
is a local theory, but not an Eliashberg theory. Indeed, in 
our case, the interaction is larger than the Fermi energy 
and the presence of the momentum cutoff in the bosonic 
propagator is crucial. This distinction becomes particu- 
larly important if we compare our result for p s with the 
conventional Eliashberg theory. There Ep is the largest 
scale in the problem, even if the "local" interaction d(ui) 
scales as 1/u or even faster (as, e.g., 1/u 2 for phonon 
superconductors). Once Ep is the largest energy scale, 
p s /kp is always larger than T pa j r , and phase fluctuations 
are weak. Indeed, according to Eq. t\'d 7|) . p s scales as 



A 



Ps 



E(w ~ A) ' 



(47) 



where, as before, p° s ~ Epkp is the stiffness of the weak 
coupling limit. The ratio A/E can be quite small if the 
pairing occurs in the quantum-critical regime and in- 
volves near-massless bosons. In particular, for phonon 
superconductors, when the Debye frequency ajp> is much 
smaller than electron-phonon interaction u, A ~ T pa j r ~ 
u (see Ref4&), and E(w ~ u) ~ u 2 /ui D » A. Then 
Ps ~ Ps( w -d/ u ) *C Ps- Still, the condition that Ep is 
the largest energy scale implies that E(w) < Ep, i.e., 
u 2 /ujb < Ep. Then, even though p s is reduced from its 
weak coupling value, it still holds that 



Ps/ki 



T, 



Epujd 



pair 



pair • 



(48) 



This implies that the exchange coupling in the corre- 
sponding XY model is still larger than the onset tem- 
perature for the pairing. As a result, within Eliashberg 
theory one can expect at most modest changes in the 
transition temperature due to phase fluctuations. In our 
case, we remind, at strong coupling Ep is no longer the 
largest energy scale in the problem, and the 1/uj form of 
d{uj) in the strong coupling limit emerges once one im- 
poses a cutoff in the integration over bosonic momenta. 



C. longitudinal gap fluctuations 

In previous subsections we discussed the role of phase 
fluctuations. They are sufficient to destroy superconduct- 
ing order between T c and T pa ; r . There also exist, how- 



ever, longitudinal fluctuations of the pairing gap, and it 
is instructive to consider how strong they are. 

Longitudinal gap fluctuations generally reflect how 
shallow the profile of the free energy with respect to de- 
viations of A(cj) from its equilibrium value is. A shallow 
profile implies that the superconducting order is weak as 
different A(u) have almost the same condensation energy. 
A situation with a shallow profile emerges when, in real 
frequencies, the attractive part of Red sc (o;) is weak, and 
the pairing predominantly comes from lmd sc (io). The 
imaginary part of a "local" interaction describes purely 
retarded interaction between fermions. This interaction 
then does not contribute to the superconducting order 
parameter, which is an equal time correlator. Accord- 
ingly, the slope of the free energy is determined only by 
a weak Ked sc (uj). 

A simple estimate of the energy scale at which longitu- 
dinal gap fluctuations become relevant can be obtained 
by analyzing the form of Red sc (u>). Converting Ea.l|41|) 
to real frequencies yields 



d sc (0) = | log M 1 



<>- 



(49) 



u O,sc) 



We see that Ked sc (uj) remains attractive up to a fre- 
quency u}q. sc /\/2, and is repulsive at larger frequen- 
cies. This means that frequencies above Wo.sc/VS do 
not contribute to the superconducting order parame- 
ter, although they do contribute to the pairing itself via 
Imd sc ((jj). This in turn implies that longitudinal gap fluc- 
tuations become strong at T > luq^ sc ~ Ep/^/g. Com- 
paring this result with Eq. Ij^6(l , we see that in our strong 
coupling limit, phase and amplitude fluctuations of the 
gap are equally important, as the corrections to the su- 
perconducting order parameter from both fluctuations 
become O(l) at T ~ T c ~ Ep/ \fg. One can equally ar- 
gue that T c <C T pa i r is the result of strong phase fluctu- 
ations, or the result of soft longitudinal gap fluctuations 
brought about by the absence of a repulsive component 
of Hed sc (u>) at fl > wq,sc- 



VI. MIGDAL PARAMETER 

As we discussed in the Introduction, the coupled equa- 
tions JSJ for the pairing vertex and fermionic and bosonic 
self-energies are valid if vertex corrections and the mo- 
mentum dependent part of the self energy can be ne- 
glected. In case of electron-phonon interaction, this ap- 
proximation was justified by Migdal 5 . Below we evaluate 
the leading corrections to our local theory, both in the 
Eliashberg limit, and at strong coupling. For definiteness, 
we focus on vertex corrections 5T of the total vertex 



(l + ST) 
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k+q, oo+£2 



iltyp relevant to the pairing problem. In what follows we 
evaluate ST for weak, intermediate and strong coupling, 
and specify in each case the relevant bosonic momentum 
and frequency. 



A. Vertex corrections at weak coupling 



k, co 



FIG. 4: Leading correction of the fermion boson vertex 
Ffc i9 (uj, fl). The solid lines stand for the fermionic propagator 
and the wiggly lines for the bosonic propagator, respectively. 



Generally, the correction to the interaction vertex be- 
tween fermions and gapless bosons, depend on the in- 
terplay between the bosonic momentum and frequency. 
In particular, Ward identities imply that vertex correc- 
tions in the limit of vanishing bosonic momentum are 
of the same order as the fermionic self-energy, and not 
necessary smalls. However, for the pairing problem, we 
need to analyze vertex corrections for typical bosonic en- 
ergies, ft typ , and for typical bosonic momenta, qt yp , that 
contribute to the pairing. 

The leading vertex correction in the normal state is 
presented in Fig I VII and is given by: 



8T q (uj,fl) = — / D k - k ,(u-u')G k ,{u)') 



(50) 



Here, uj is the external fermionic frequency. In principle, 
ST depends on two momenta - the bosonic momentum 
q and the external fermionic momentum k. However, 
the dependence on k is weak and thus irrelevant, and we 
will neglect it. Performing the momentum integration in 
(I5t)ll . we obtain 



ST q (uj, ft) 



3g J Q d (uj + lu') doj' 
yJw + {Nv F q) 2 



(51) 



The factor TV in the denominator is a consequence of the 
rescaling that we performed in Sec. II. In the limit q — > 



5T g ^ (uj,fl) = JL I d{uj')duj' (52) 

£ (ft + uj) - £ (uj) 
^ift ' 

This is the Ward identity relating the homogeneous ver- 
tex with the self energy. We see from (|52^l that static 
vertex corrections do not depend on N and are not small 
at moderate and strong coupling. The situation, how- 
ever, changes when we evaluate ST q (uj, ft) at <7t yp and 



We first consider the limit of weak coupling, g <C 
1. The typical bosonic energy ftt yp is of order T c ~ 
(Ep/g)e^ 7T ^ 2 ^K On the other hand, typical momenta 
<7ty P are obtained from the condition that the momentum 
and frequency dependent term in the bosonic propagator 

1/3 

D(q, ft) are of the same order, i.e. qt yv — (7Slt yp ) ■ 
Together with Ea. lfTUI) for 7 this yields 



vpgtyp ^ E F e 



-w/(6v?) 



Comparing Q and v F q, we see that 



VFqtyp 



n 



typ ge 



,t/(3Vs) 



> ft 



typ- 



We then obtain from Ea. (|51|l 



5T 



Nv F q tTP 



^g 



(53) 



(54) 



(55) 



We see that vertex correction is exponentially small for 
small g. The extension to large N is in fact not needed 
as vertex corrections are already negligible. 



B. Vertex corrections at intermediate coupling 

At intermediate g = 0(1), ujq and Ep become of the 
same order, i.e., there is only one characteristic energy 
scale in the problem. Vertex corrections are 0(1) for 
N = 1, but are still small in 1/N if we extend the theory 
to large N. 



C. Vertex corrections at strong coupling 

A naive expectation would be that vertex corrections 
gradually increase with g and eventually overcome the 
overall smallness in 1 /N. This would invalidate our local 
theory at sufficiently large g. It turns out, however, that 
vertex corrections freeze at 0(1/ N) and do not grow with 
g. The saturation originates from the form of the bosonic 
propagator D(q, ft), which is determined by the self en- 
ergy ~ 7 g • As 7 by itself scales with the boson- fermion 
coupling, D(q,uj) scales inversely with g and cancels out 
the overall factor g in Eal5D 

To see this explicitly we note that at strong coupling, 
qtyp is of order kp, hence vpqt yp is of order Ep. Typical 
frequencies ft tyv are also 0(Ep). The external fermionic 
frequency w in Eq. I|51|l is also of order of the Fermi 
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energy. Evaluating the vertex correction diagram using 
d (Q) ~ jjpy , we obtain for these q typ and fi t yp 



it follows from our analysis that the pairing sets in at 



pair 



ST 



r 1 
9 Jo 



Ep-\-u) 



■dw' 



1 

N 



(56) 



We see that vertex corrections indeed do not depend on 
g and remain small (at large N) for arbitrary strong cou- 
pling. 

At large N, vertex corrections remain small for all val- 
ues of g, i.e., our local theory is valid both in the weak 
and the strong coupling limit. 



VII. OTHER SYSTEMS 

As we discussed in the introduction, the problem dis- 
cussed in this paper describes p— wave superconductiv- 
ity in condensed matter systems close to a ferromagnetic 
quantum critical point with Ising symmetry™, and super- 
conductivity due to color magnetic interactions in dense 
quark matte»22i2L2&. As pointed out above, for color su- 
perconductivity, antiparticlc pairing that has been ne- 
glected in Eqs. JSJ), J7J), may come into play at strong 
coupling 41 . We nevertheless believe that the main re- 
sults of this paper are still relevant to this case. 

However, the results obtained in the strong coupling 
limit are much more general. The key aspect of our the- 
ory is that the boson propagator -D(qii, q±, f2) becomes 
completely local above a certain frequency. We consid- 
ered a particular case when boson dynamics is set by 
Landau damping, and 



(57) 



where luq is a characteristic upper cut off scale of the 
boson system. However, our results will be equally valid 
for any D(O) in the form D(ft) oc 1/ |0| 7 with 7 > 1. 

One of the possible applications of our strong cou- 
pling result is pairing by antifereromagnetic fluctuations, 
which has been discussed in great detail in the context of 
cuprate superconductors^. That analysis was, however, 
performed within the low-energy spin-fermion model, 
which is only valid at u < Ep, i.e., in the weak coupling 
limit, as we defined it in this paper. In the cuprates, 
as is well known, the Hubbard interaction is comparable 
to Ep, i.e., g > 1. Our finding that a large pseudogap 
regime with phase incoherent pairs is inevitable at strong 
coupling is quite intriguing in view of the numerous ob- 
servations of pseudogap physics in this class of materi- 
als. To make a connection to other studies of pseudogap 
within Hubbard model, we note that our effective inter- 
action u ~ U 2 /t and Ep ~ t, where U and t are the local 
Coulomb repulsion and tight binding hopping element, 

respectively. Then g = (-j) and u>o ~ Ep/g oc jj J, and 

wo, s ~ Ep/y/g ~ J, where J = ^j- is the antiferromag- 
netic exchange interaction between spins. Accordingly, 



t, and the pairing gap is A ~ t, while coherent 
superconductivity occurs at T c oc J. We note in this 
regard that numerical investigations of variational wave 
functions designed to cover the strong coupling limit 57 
do indeed yield a zero temperature gap A ~ t and a 
considerably reduced superfluid stiffness ^.SiSLSI. Also, 
our result T c ~ J is consistent with the observation that 
T c in underdoped cuprates scales with the neutron peak 
frequency. The latter turns out to be of order J, if, e.g., 
one extends the exciton scenario for the resonance peak 
to the strong coupling limit 61 . 

There exists some similarity between our results and 
those obtained for the crossover from BCS-type behavior 
at weak coupling and Bose Einstein condensation (BEC) 
of pairs at strong coupling 62,63 - 64,65 . In particular, for 
large g, when A is or order Ep, our pair coherence length 
£ ~ vp/A becomes of the order of the typical distance 
between fermions (~ kp 1 )- This is similar to the find- 
ings for the BCS-BEC crossover 62 i 63 i 64 i 65 . An important 
distinction between the two theories is that in BCS-BEC 
crossover, the pairing interaction is static, while in our 
case it is dynamic and strongly retarded. The transition 
at large g in our case should therefore not be considered 
as condensation of almost free bosons. 



VIII. CONCLUSION 

In this paper we considered pairing of 3D fermions, due 
to an exchange of massless bosons. The model we con- 
sidered describes p— wave superconductivity in itinerant 
fcrmionic systems close to a ferromagnetic quantum crit- 
ical point with Ising symmetry, and superconductivity 
due to color magnetic interactions in quark matter. 

At weak coupling, we find that an exchange of a mass- 
less boson enhances T c compared to the BCS expectation 
T^ cs tx e~9, and the actual transition temperature is 
T c ~ ujQe 2 v^, where loq = Ep/g is a characteristic en- 
ergy scale of the bosons. This result agrees with previous 
calculations 32 ! 37 ! 38 . 

At strong coupling, we find that pairing emerges at a 
temperature T pa ir ~ 0.06-Ef, which is only numerically 
smaller than the Fermi energy. In addition, we find that 
the phase stiffness behaves as p s ~ Tp^/^Jg, i.e. the 
typical energy scale where phase fluctuation become im- 
portant is parametrically smaller than T pa h-. The small 
phase stiffness implies that coherent superconductivity 
only emerges at T c ~ T-p^/y/g <§; T pa i r . In between T c 
and Tp a i r the system displays pseudogap behavior with 
preformed pairs. The width of the pseudogap regime 
widens as g grows. 

We further argued that several aspects of the strong 
coupling solution, particularly the existence of two tem- 
peratures Tp a i r and T c , are valid for a much larger class of 
problems in which boson propagator becomes completely 
local above a certain frequency. The existence of a cut 
off energy scale, which corresponds to the finite lattice 
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constant in the condensed matter context or to the large 
density in case of color superconductivity, is crucial for 
the relevance of phase fluctuations. For the cuprates, our 
results imply that T pa i r and the pairing gap A scale with 
the hopping integral t, while T c scales with the exchange 
interaction J. 



it must hold that y (b) = for any 7^1. For real b the 
condition y (b) = cannot be fulfilled. However, for a 
complex b = a + i/3, we find that the imaginary part of 
y (b) vanishes if a — ^, i.e. if y Q + i/3) is purely real. 
Using this fact and substituting b = 1/2 + i/3 into l|A5() . 
we obtain that /3 is determined from 
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This equation is easily solved graphically and yields 
±0.7923(2). Therefore 



hy (uj) = A lu-?-^ + A*uj-^ +if) . 



(A7) 



The overall constant A is chosen such that h 7 (uj) is real. 

In order to show that Ea. (|A7Jl is indeed the solution 
of Ea. H32|) we discuss more carefully why lim 7 _>i /i 7 (u>) 
gives the desired $ (uj). Eq. ll3*2l can be re-expressed as 



lim 

7—1 



1 



7 



E p 



duj'<& (uj 1 ) fc 7 (u),u)') 



yi-7 



CO, 



1-7 



(A8) 



APPENDIX A: PAIRING WITH 



KERNEL 



In this appendix we obtain the solution of the Eliash- 
berg equation, Ea. (|3*2")l . in the strong coupling regime. 
The analysis of a pairing problem with a kernel 



1 rr -l>j 
\u—oj I 

nontrivial and it is useful to solve a more general problem 
first 66 «. We consider: 



This equation coincides with IjAlfl if we neglect Wg -7 in 
the denominator. We now recall that at strong coupling, 
Tpair > wo- Then for all u>' > T in (|X8|) holds that J > 
u>o, and we can safely neglect luq ~ 7 for any 7 7^ 1. 



APPENDIX B: LINEARIZED GAP EQUATION 
FOR DISCRETE MATSUBARA FREQUENCIES 



w .!^r*5^M (ai) 



(u-)'-i 



with 



fey (W, Uj') 



(A2) 



and argue below that $ (ui) ~ lim 7 _>i /i 7 (to). 

The integral equation (|A1|) is scale invariant suggesting 
a power-law solution 

h-y (ui) = Auj- h . (A3) 

Inserting this ansatz into Eq. (|Alj) we obtain 

I -- [ s ,1/1 1 
' 1 dt- 



1 



t i+i-7 vi* l r it +11 



(A4) 



where t = — . This determines the exponent 6(7). The 
integral over t can be performed explicitly. In the limit 
where 7 is close to 1, Eq. (|A4(I reduces to 



1= 1 + (1-7)2/ W 



(A5) 



where y (b) — je + i> (b) — f tan (^f\ , and i(j (b) is the 
di-gamma function. While b is undetermined for 7 = 1 



In the computation of T pa j r in the main text we im- 
posed the lower cutoff in the zero-temperature equation 
for the pairing vertex. In the weak coupling limit, this 
procedure was shown earlier— to yield the same T pa i r 
(modulo a numerical prefactor), as one would obtain by 
performing an explicit summation over discrete Matsub- 
ara frequencies. In this appendix we demonstrate that 
the same is true in the strong coupling limit. 

The most straightforward way to analyze the linearized 
pairing problem is by considering the gap function 



An = 



Z (ui n ) ' 



(Bl) 



where Z (u> n ) = 1 - 



and determine the temperature 



at which A„ 7^ for the first time. The advantage of ana- 
lyzing A instead of $ is that the corresponding linearized 
gap equation does not explicitly contain the fermionic 
self-energy and can more easily be solved numerically. 
The equation for A„ is straightforwardly obtained from 
the equations for $ (w„) and S (ui n ): 



An = 3ngT V ( ^ - ^ 



xsig 



;n (u) m ) d (uj m - 



•n 



(B2) 
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In the limit where T pa ; r > ojq, relevant \u> n \ > ojo, we can 
approximate d(fl n ) by d(Q n ) ~ | . Then 

» _ -Ef v-* / A m A„ \ 

" ~ 4^T pair ^ V 7i7+l/2 ~ 71 + 1/2; 

^ 7717=71 

x Slgi ; (m+ ; /2) . (B3) 

1 771 — 7l| 

The summation over discrete Matsubara frequencies is 



convergent for large n. Thus, no regularization or cut off 
is needed to solve for A m . The only dimensionless pa- 
rameter in (|B.3J| is Ep/T pa x T , and the non-trivial solution 
of 1B3|1 . if it exists, appears at T pa j r ~ Ep as we found 
earlier. We verified numerically that the solution of Eq. 
(|B3j) does indeed exist at T pa i r ~ Q.06AEp. This is even 
quantitatively close to the estimate T pa ; r ~ 0.0676-Ef ob- 
tained in the main text. 
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